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Abstract. Generation of magnetic field energy, without mean field generation, 
is studied. Isotropic mirror-symmetric turbulence of a conducting fluid amplifies the 
energy of small-scale magnetic perturbations if the magnetic Reynolds number is 
high, and the dimensionality of space d satisfies 2.103 < d < 8.765. The result does 
not depend on the model of turbulence, incompressibility and isotropy being the 
only requirements. 

91.25. C. 



Solar and galactic magnetic fields are generated by turbulent motions of their 
electrically conducting quasifluid constituents. Generation of large-scale or mean 
magnetic fields is described by the kinematic dynamo theory of Steenbeck, Krause 
and Radler [1]. Small-scale magnetic fields can come from the large-scale fields 
because of turbulent stretching and twisting and play an important role in satura- 
tion of the dynamo generated mean fields [2-4]. It was noted long ago [5-7] that 
small-scale fields (SSF) can also be generated in the absence of mean fields. For 
example, mirror-symmetric turbulence is incapable of the mean field dynamo, but 
can amplify SSF perturbations. A nice qualitative picture of this process was given 
by Moffatt [8] . When the magnetic Reynolds number is large, the magnetic field is 
frozen into the fluid. The field magnitude grows as the infinitesimal line element, 
that is exponentially. However, the characteristic wavenumber of the magnetic per- 
turbation also grows exponentially and finally reaches the dissipation scale. It is a 
delicate question which one of these effects wins on the average. 

In many astrophysically interesting cases the viscosity is much larger than the 
resistivity, and the turbulent velocities cut off at a scale much larger than the 
scale at which resistivity destroys magnetic field. The SSF dynamo exists on scales 
between the viscous and resistive cut offs. In high Reynolds number turbulence the 
SSF dynamo growth rate exceeds the growth rate of the mean field dynamo by a 
factor of Re 1 / 2 . Thus, if SSF dynamo does operate it should be important for the 
large-scale dynamo saturation [3]. Also of interest are the consequences of the SSF 
dynamo for the applicability of the kinematic mean field dynamo as discussed by 
Kulsrud and Anderson [9]. In this letter, we show that the SSF dynamo indeed 
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operates for spatial dimension d in the interval 2.103 < d < 8.765. The result is 
model-independent, in particular turbulence is not assumed to be 5-correlated in 
time. 

To get this result, we use a zero-dimensional representation of the induction 
equation. Consider a small element of fluid with dimensions much smaller than the 
smallest scale of the turbulent flow, but much larger than the scale of the magnetic 
field. We may then approximate the velocity field of the fluid element by a Taylor 
expansion around an arbitrary center r . Thus, v a = fo a + V^fe(r& — r &) + .... Here 
V is the rate of strain tensor moving with the fluid element, V a b = d a Vb- Note that 
the rate of strain tensor is dominated by the smallest scales close to the viscous cut 
off [9]. We assume incompressibility, trV = 0. Then the magnetic field evolution is 
given by equations for the field B and the wavenumber k [10]: 



here V* is the transpose of V; the molecular magnetic diffusivity will be taken into 
account later. Equations (1),(2) are applicable for SSF dynamics only, at large 
length scales they are invalid. More precisely, (1), (2) are always exact, being just 
equations for co- and contravectors. However, if the length scale of the magnetic 
field becomes comparable to the length scale of the turbulent velocity field, it takes 
infinite number of such equations to describe the time evolution of the magnetic field. 
This infinite system of equations is equivalent to the standard induction equation 
dtB = V x (v x B). Formal solution to this linear equation can be written, but we 



B = VB, 



(1) 



k = -v% 



(2) 
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do not know how to average the resulting magnetic energy. The linear (kinematic) 
dynamo is a non-trivial problem. 

The system (1),(2) is much simpler than the standard induction equation, all 
the relevant properties of turbulence being concentrated in the stretching tensor 
V. We wish to know whether the energy of the SSF grows or decays for a random 
stretching rate V with given statistics. 

As an introduction to the problem, consider the Gaussian 5-correlated turbu- 
lence. Then turbulence is specified by the pair correlator of the stretching tensor. 
From isotropy and incompressibility, 

< V ab V cd >= 2~f8(t)[(d+ l)8 ac 8 bd - 5 ab S cd - 8 ad 8 bc ), (3) 

where 7 is a constant with dimension of frequency. The physical meaning of 7 is the 
characteristic stretching rate, 7 is determined by the smallest eddies. Let P(k, B) 
be the probability distribution for the pairs (fc, B). Since V of (1),(2) is 5-correlated 
in time, one can derive a Fokker-Planck equation for P using standard procedures 
(see eg. [11]): 

d t P = (l/2)[d ka d kb D k a k b + 2d ka d Bh D k £ + d Ba d Bb D* b B ]P, (4) 
where the mean velocities vanish and the diffusivities are 

D$ = >y[(d+l)k 2 6 ab -2k a k b ], (5) 

Da* = l[(d+l)B 2 8 ab - 2B a B b ], (6) 

D k J? = ~l[(d + l)k b B a - k a B b - k c B c 8 ab }. (7) 
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For our purposes, it suffices to now the magnetic energy spectrum W(k) denned as 
a S 2 -moment of the probability distribution P. From (4-7), taking into account 
solenoidality k c B c = 0, and assuming isotropy, one gets the evolution equation for 
the magnetic energy spectrum 

W = (d- l)k 2 W" + (d 2 - 5)kW + 2(d 2 -d- 2)W, (8) 

where prime stays for the modulus k derivative, and time was normalized to exclude 
a factor proportional to 7. For d = 3 this equation coincides with those obtained 
in [7,9]. 

Boundary conditions for (8) is the immediate question. First of all, molecular 
magnetic diffusivity will eat up all the modes with very large k. Thus, the boundary 
condition at k = 00 is simply W = 0. It is assumed that the limit k max — > 00 is 
taken after all other limits. The order of performing limits can be important. For 
example, magnetic energy E = J dkW in two dimensions without resistivity grows 
exponentially, 

lim lim E = 00. 

At the same time, magnetic energy finally decays in two dimensions if the resistivity 
is positive, thus 

lim lim E = 0. 

This follows from the Zeldovich antidynamo theorem, and our calculations will 
reproduce this result. 

Secondly, the very model (1),(2) is not applicable at small k, where the turbu- 
lence is not reducible to the corresponding rate of stretching. In fact, at large scales 
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all we have is a usual turbulent magnetic diffusivity. The magnetic field is thrown 
out of the large-scale region into turbulent scales. The corresponding boundary 
condition should be W = at some k = k m i n . From what follows and because 
k m i n is small, this boundary condition is equivalent to W = at k = 0. 

The eigenmode of (8) is simply k v , and the growth rate is 

7 (y) = (d- l)y(y - 1) + (d 2 - 5)y + 2(d 2 -d-2). (9) 

To satisfy the boundary conditions, one has to have a pair of solutions with the 
exponent y = Rey ± iO giving the same growth rate 7 (then a linear combination 
of these modes will vanish at logk = ±00). This means that y is determined by the 
requirement 

i'(y) = 0. (10) 

Calculating y and then calculating the corresponding growth rate, we get the SSF 
dynamo (7 > 0) for dimensionalities d satisfying 

d{d- l)(9-d) > 16, (11) 

or 2.103 < d < 8.765 as advertised. 

Certainly we are mostly interested in d = 3 case, which happens to be a dy- 
namo case according to (11). However, working in d dimensions is useful for what 
is coming. Suppose, that the ^-correlation supposition is dropped. Then the coef- 
ficients in the diffusion equation (8), and in fact the very form of the equation will 
change. Perhaps we will still have SSF dynamo in some dimensionalities interval, 
but the lower critical dimension might turn out to be greater than 3, and there will 
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be no SSF dynamo in three dimensions. As Vainshtein [12] puts it, the SSF dynamo 
is a quantitative rather than a qualitative problem. It was our original idea to show 
that the critical dimensions do depend on the correlation properties of turbulence, 
and the 5-correlated result of [7,9] is not reliable. To our surprise, we found that 
the critical dimensions are structurally stable and probably even universal as far 
as the turbulence is incompressible and isotropic. We give a short account of our 
investigations in what follows. 

The idea is to solve (1),(2) for a long rather than a short time interval. The 
solution is 

B = UB , (12) 

k = (U t )~ 1 ko. (13) 

Here U is the evolution matrix given by U = VU ', U(t = 0) = J; in terms of Euler r 
and Lagrange a coordinates U = dr/da. If time of integration is long enough, the 
matrix U can be written as R1DR2, where R are random rotation matrices, and 
D is diagonal. Physically speaking, after several correlation times, the orientation 
of the fluid element is randomized. Now suppose that at time t=0 the magnetic 
energy spectrum was isotropic, Wo = Wo(k), where k now means the magnitude. 
Using randomness of the rotation matrices R, one can calculate the new spectrum 
W{k) after the transformation (12), (13) as 

W(k) =< iDieWD-^-V+VWodDklk) >, (14) 
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where < ... > means averaging over the angles of the perpendicular <i-dimensional 
unit vectors k and b. The integral transformation (14) generalizes the Fokker-Planck 
equation (8). Equation (14) is the main result, it can be derived as follows. 

If Pq was the probability distribution of pairs (k, B) before the transformation 
(12), (13), then the new pdf P after the transformation is 

P(k,B) =P (U%U~ 1 B). (15) 

Because of isotropy, the magnetic energy spectrum W(k) can be calculated as 

W(k) = J B d+1 dB < P(kk,Bb) >, (16) 

where < ... > is the angle average. Now plug (15) into (16) 

W(k) = J B d+1 dB < P^kU^k.BU- 1 })) > . (17) 

As a consequence of isotropy and solenoidality, the probability distribution has a 
form 

P (k,B) = Q(k,B)6(k-b), (18) 

and (17) can be written as 

W(k) = J B d+1 dB < Q{\U l k\k, |t/" 1 6|S)5(A; • S/flE^fcHEr^D) > . (19) 

Now change the integration variable in (19), B — > B/\U~ b\ and get 

W{k) = j B d+1 dB <\U t k\\U- 1 b\-( d+ ^Q(\U t k\k,B)5(k-b) > . (20) 

Recalling (18) and the definition of Wo, we write (20) in a form (14). Only the 
diagonal part of U matters because of the angle averaging. 
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A nice property of (14) is that the eigenmode is still k y . The amplification 
factor after a single transformation (14) is given by 

F(y) =< | J Dfc|^+ 1 )| J D- 1 6|-( d+1 ) > . (21) 

Write the diagonal part as D = exp diag(Xi/2, Ad/2), the incompressibility im- 
plies Ai + ... + Ad = 0. Then, 

\Dk\ 2 = J2 eXik j, ( 22 ) 

3 

iD-i&l 2 = J2 e ~ Xjb2 3- (23) 

3 

Like before, to satisfy the boundary conditions one chooses the exponent y in such 
a way as to minimize V: 

T'(y) = 0. (24) 

Thus, the necessary and sufficient condition for the SSF dynamo is that V be greater 
than 1 for arbitrary y and Xj. Since V is equal to 1 at Ai = ... = Ad = 0, this point 
should be a minimum for the dynamo case. Calculating (21) in the vicinity of zero, 
one gets 

F(y) = l + [ 4 rf(rf-l)(rf+2)]- 1 (^A J 2 )[(rf-l)y( 2/ -l) + (rf 2 -5) ?/ + 2(rf 2 -rf-2)]. (25) 

3 

Comparing to (9), we see that V is greater than 1 in the vicinity of of the point 
A = if (11) is satisfied. Thus the structural stability of the result (11) is proved. 
We suspect that (11) ensures that V is greater than 1 everywhere, not only in the 
vicinity of the zero point. We do not have a proof. However, numerical calculations 
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of the integral (21) in three dimensions confirm that A = is a global minimum 
and T is greater than 1 everywhere else. 

To conclude, we have shown that SSF dynamo operates in three dimensional 
isotropic turbulence. The result seems to be model independent and is definitely 
structurally stable. The SSF dynamo is of great importance to the understanding of 
the growth of large scale fields. Galaxies and upper layers of stars are characterized 
by R m 3> Re which is the necessary condition for a SSF dynamo to be present. The 
large growth rate of the SSF dynamo - it is Re 1 / 2 bigger than the mean field growth 
rate - ensures that self consistency effects first enter at the small scales, when the 
large scale fields are negligible [9] . Clearly large scale fields are observed in galaxies 
and stars, and the dynamo process must involve considerable evolution where the 
small scale fields are dynamically important but the large scales are not. The SSF 
dynamo saturation mechanisms and the effects of SSF dynamo on the conventional 
mean field dynamo should be subjects of further studies [9,13]. 
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